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transmission resonance that can persist when the broadening of individual meta-atom resonance fre- 
quencies is roughly one half the frequency characterizing the split ring asymmetry. We also find that 
larger degrees of inhomogeneous broadening can drastically alter the cooperative response, destroy- 
ing this resonance. The reduced effect of cooperative response due to inhomogeneous broadening 
may find applications in producing metamaterial samples that more closely mimic homogeneous 
magneto-dielectric medium with well-defined susceptibility and permittivity. 
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I. INTRODUCTION 

There exists increasing experimental evidence that 
metamaterial samples, consisting of arrays of sub- 
wavelength scale nano-structured circuit elements, can 
be prepared in the limit where collective interactions be- 
tween the resonators play a vital role in determining their 
electromagnetic (EM) responses. For example, experi- 
ments on two-dimensional (2D) arrays of closely spaced 
asymmetric split ring (ASR) meta-molecules have indi- 
cated the presence of a high-quality transmission res- 
onance, with a dramatic sensitivity of the resonance 
linewidth to the number of ASRs in the sampled The 
transmission resonance was also observed to depend 
strongly on the spatial distribution of the mutually in- 
teracting resonators. Where a sufficiently large, reg- 
ular array of ASRs exhibits a high quality transmis- 
sion resonance* 1 - introducing disorder in the elements' 
positions destroys the observed spectral transmission 
window j2i£ further indicating the role collective EM inter- 
actions can play in metamaterial dynamics. We recently 
showed 4 ^ 5 - that closely spaced ASR metamolecules inter- 
acting via a resonant EM field exhibit collective eigen- 
modes with strongly suppressed resonance linewidths. 
The cooperative response yields the characteristic fea- 
ture in the experimentally observed enhanced quality fac- 
tor of the transmission resonance in Ref. [l|. Numerically 
analyzing the properties of collective modes with narrow 
radiative resonance linewidths provided a physical expla- 
nation of this phenomenon with an excellent agreement 
between the simulations and the measurements. 

In our previous study 5 - showing how the transmission 
resonance observed by Fcdotov et al* depends on the 
linewidth of a particular collective mode, all unit-cell res- 
onators were assumed to respond to EM fields identically. 
In the preparation of metamaterial samples, fabrication 
defects, however, may in general lead to variation in the 
geometry of individual resonators. The current oscilla- 
tions supported by the unit-cell resonators would there- 
fore possess slightly different resonance frequencies, re- 
sulting in inhomogeneously broadened metamaterial ar- 



rays. Inhomogeneous broadening changes the conditions 
of the resonant interaction processes. This may impair 
the collective, coherent phenomena which are potentially 
important in several applications and physical effects of 
metamaterials such as lasing^ and providing precise con- 
trol and manipulation of EM fields on a sub- wavelength 
scale, as theoretically proposed in Ref. and experimen- 
tally observed in Ref. @ (for other related studies, see for 
example Ref. d) . 

In this work, we study how uncontrolled inhomoge- 
neous broadening of plasmonic resonators limits the ob- 
servation of collective phenomena in metamaterials. We 
will examine under what conditions these collective ef- 
fects can still persist and potential implications of inho- 
mogeneous broadening on metamaterial applications. In 
particular, we consider a rectangular 2D array of ASR 
meta-molecules whose resonance frequencies are stochas- 
tically varied. We numerically evaluate ensemble av- 
erages of the EM response of the metamaterial over 
the stochastic distributions of the resonator properties. 
These show how increased inhomogeneous broadening 
inhibits the response of the coherent collective modes 
responsible for the transmission resonance observed in 
Ref. [l|. Moreover, we find the effects of interactions be- 
tween different discrete resonator elements, which result 
in the narrowing of the radiative resonance linewidth, are 
diminished as a function of an increasing inhomogeneous 
broadening. Our results therefore illustrate how main- 
taining uniformity in the fabrication process is essential 
in designing new metamaterial based devices whose ap- 
plications rely on strong interactions between the res- 
onator elements and on a cooperative response. 

The multiple scale spatial structure associated with 
nano-fabricatcd resonators in metamaterial arrays, along 
with the wave nature of scattered EM fields, poses a the- 
oretical challenge for studies of the response of these sys- 
tems to resonant EM fields. Interactions resulting from 
recurrent scattering events, in which a field is scattered 
more than once by the same resonator, frequently play a 
crucial role in the cooperative system responses! 4 ' 5 ' 10 " — 
While in an infinite, regular lattice, the translational 
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symmetry can be exploited to calculate approximate lo- 
cal field corrections in a medium of discrete scatterers^, 
recurrent scattering processes are generally more difficult 
to model in finite-sized systems with complex geometries. 

However, since inhomogeneous broadening reduces the 
cooperative effects arising from recurrent scattering, en- 
gineering a controlled amount of inhomogeneous broad- 
ening into the metamaterial, e.g. via geometrical vari- 
ation of the resonators, may provide a practical means 
to produce samples that are easier to design and theo- 
retically analyze. In particular, we find that with an in- 
creasing inhomogeneous broadening the response of the 
system approaches that of standard continuous medium 
electrodynamics. Reducing cooperative effects is po- 
tentially important because several metamaterial appli- 
cations, such as diffraction-free lenses formed from a 
medium with a negative refractive index,— ~— are sim- 
plest to realize with a well-defined electric susceptibility 
and magnetic permeability that, in many systems with 
complex geometries, are only approximately achieved^ 

The effects of inhomogeneous broadening have previ- 
ously been examined, e.g., by Gorkunov et al on the 
bulk properties of left-handed materials in periodic in- 
finite lattices^ and have been experimentally observed 
by Gollub et al*£- In our study of ASR resonators we 
evaluate collective modes of a finite lattice in order to 
investigate the effects of the inhomogeneous broadening 
on the experimentally observed sample size dependent 
transmission resonant lincwidth narrowingji 

In our analysis we employ a general theoretical for- 
malism of collective interactions between a discrete set 
of plasmonic resonators, or meta-atoms, mediated by the 
EM field that we developed in Ref. 0. In the model, we 
assume each meta-atom exhibits a single mode of cur- 
rent oscillation that possesses appropriate electric and 
magnetic dipole moments. Each meta-atom responds to 
EM fields exhibiting a specific resonance frequency and 
coupling strength that are determined by its character- 
istic design. Starting from the Lagrangian and Hamil- 
tonian formalism describing the interaction of the EM 
field with polarization and magnetization densities cre- 
ated by a charge distribution, we then derived the cou- 
pled dynamics of the EM fields and the meta-atom dy- 
namic variables^ In a collection of meta-atoms, inter- 
actions with the EM field mediate a dynamic coupling 
between the meta-atoms and determine the collective dy- 
namics within the ensemble, resulting in distinct collec- 
tive modes with corresponding resonance frequencies and 
lincwidths. The analysis of collective response in terms 
of discrete resonators also points to the direction of an 
interesting analogy between resonators and a system of a 
cloud of atoms. In atoms the electron transitions driven 
by an EM field create an electric dipole moment, while in 
the case of circuit elements the oscillating current gener- 
ates both the electric and magnetic dipole moments. The 
model of Ref. has previously been successful in provid- 
ing an excellent agreement between the theory and ex- 
perimental observations of cooperative transmission res- 



onance linewidth narrowing of ASR metamolccules^ 

The remainder of this article is organized as follows. 
We summarize the theoretical formalism we employ to 
describe collective interactions within the metamaterial^ 
in Sec.|TTJ In Sec. lIIIl we describe the fundamental build- 
ing block of our metamaterical, the ASR, in the context 
of this formalism. The main results of the article, de- 
scribing the effects of inhomogeneous broadening on the 
collective response of the metamaterial are presented in 
Sec. IIV1 and conclusions follow in Sec. [V] 

II. A MODEL FOR THE COLLECTIVE 
RESPONSE IN METAMATERIALS 

In order to incorporate the effects of strongly hetero- 
geneous metamaterial we describe the sub-wavelength 
structures of the medium as discrete scatterers^ Each 
unit-cell element, a meta-molecule, may also consist of 
sub-elements, which we call meta-atoms. While the gen- 
eral formalism of Ref. [3 allows for multipole-field ra- 
diation of the resonator unit elements, as a first ap- 
proximation here we consider each sub- wavelength-sized 
meta-atom simply as a radiating dipole and ignore its 
multipole-field contribution. Following our treatment in 
Ref. |j, we assume that each meta-atom j, with its posi- 
tion vector defined by Yj , supports a single cigenmode of 
current oscillation. The dynamics of this current oscilla- 
tion are determined by the dynamic variable Qj(t) with 
units of charge. Each meta-atom exhibits an electric and 
magnetic dipole moment. These may be expressed as 

dj = Qjhjdj , (la) 
'" = IjAj&j , (fb) 

respectively. Here Ij(t) = dQj/dt denotes the current, 
and the directions of the dipole moments are specified 
by the unit vectors dj and rhj with proportionality co- 
efficients hj and Aj (with units of length and area, re- 
spectively) that depend on the specific geometry of the 
resonators. We assume the meta-atoms are designed 
such that the electric quadrupole and higher order mul- 
tipolc contributions to the meta-atom dynamics can be 
ignored. Although each meta-atom possesses only elec- 
tric and magnetic dipoles, a meta-molecule of two or more 
meta-atoms in our model would exhibit a non-vanishing 
quadrupole field. While, in general, this quadrupole con- 
tribution is inaccurately represented in the dipole ap- 
proximation, in the case of the ASR meta-molecules con- 
sidered in Rcfs. Q] and and in the present study, the 
generated quadrupole field is notably suppressed when 
compared to the corresponding dipolar field. This has 
been indicated by finite element simulations of Maxwell's 
equations within a single meta-molecule^ Additionally, 
the fact that the size of the meta-atoms is often com- 
parable to the spacing between them could result in a 
correction to the coupling strength between neighbor- 
ing elements obtained in the point dipole approxima- 
tion. Nonetheless, this model in the dipole approxima- 
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tion was employed in Rcf. ID to characterize the coopera- 
tive linewidth narrowing responsible for the enhancement 
of quality factor with system size observed in Ref. [l|, 
yielding excellent agreement with experimental results. 
In this section, we describe the key features of our theo- 
retical formalism that are required to describe the collec- 
tive response of an inhomogeneously broadened sample 
of ASRs to the EM field. Details of the derivation are 
presented in Ref. 0. 

We write the polarization and magnetization densities 
as a sum of their contributions from the individual meta- 
atoms 



P(r) = 

3 

M(r)=^M 3 (r), 



(2a) 
(2b) 



where the polarization and the magnetization of the res- 
onator j in the dipolc approximation read 



P,-(r,t)«d^(r-r 
M,(r,t)«m^(r- 
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(3a) 
(3b) 



respectively. 

An external beam with electric field E; n (r, t) and mag- 
netic held H in (r, t) with frequency ri impinges on the en- 
semble of meta-atoms. The incident EM field drives the 
meta-atoms, generating dipolc radiation from the oscil- 
lating electric and magnetic dipoles. The total radiation 
from the metamatcrial array is the sum of the scattered 
electric and magnetic fields from all the meta-atoms 



Eg(r,t) = ^E Sij (r,t). 

3 

H s (r,t)=^H Sj (r,t). 



(4a) 
(4b) 



where Esj(r, t) and Hgj(r,i) denote the electric and 
magnetic field emitted by the meta-atom j. The Fourier 
components of the scattered fields have the familiar ex- 
pressions of electric and magnetic dipole radiation^ 



P r 



dV 



G(r-r',n)-P+(r',fi) 



-G x (r-r',0)-M+(r',ft) 



H+ (r,Q) = - /d 3 ~' 



if 
4tt 



G(r - r', fl) ■ Mt(r', Cl) 



-cG x (r-r',f2)-Pt(r',0) 



(5) 



(6) 



where we have defined the positive and negative fre- 
quency components of a real time varying quantity V(t) 
such that for a Fourier component of frequency f2, 
V ± (Cl) = 0(±n)V(fi), and hence V(t) = V+{t) + V~(t) 
with V~(t) = [V + (f)]*. Here G denotes the radiation 
kernel representing the electric (magnetic) field emitted 



from an electric (magnetic) dipole^ The explicit expres- 
sion for the corresponding radiated field from a dipole v 
reads 



G(r, CI) ■ v = (f X v) X r- 



fa 



(fcr) 3 (fcr) 5 



[3f (f • v) — v] 

-in 



(5(fcr)v, (7) 



where f = r/r and k = fi/c. Similarly, G x (r, CI) repre- 
sents the radiation kernel for the magnetic (electric) field 
of an electric (magnetic) dipole source^ Specifically, the 
corresponding radiated field from a dipole v yields 



ikr / i 

G x (r,n).v=— 1- — 

kr \ ikr 



r x v . 



(8) 



The polarization and magnetization densities appear- 
ing in Eqs. ([5]) and ^ are themselves driven by the in- 
cident and scattered EM fields. This driving, combined 
with the scattered EM fields, yield a coupled set of equa- 
tions for the resonators and EM fields that we derived 
in Ref. 0. Current excitations in each meta-atom j pro- 
duces a field that interacts with the current and charge 
oscillations that generated it. Due to these self-generated 
fields, a meta-atom j exhibits behavior similar to that of 
an LC circuit with resonance frequency^ 



f 



(9) 



where Cj is an effective self-capacitance and the effective 
self-inductance is Lj. In this work, we consider an inho- 
mogeneously broadened sample of N ASRs. Each ASR, 
I, (I = 1 . . . N) consists of two meta-atoms whose reso- 
nance frequencies Uj (j = 21 — 1,21) are centered around 
the random frequency ujq + Xi , where Xi are independent 
identically distributed random variables. We assume the 
meta-atom resonance frequencies occupy a narrow band 
about wo, i.e. \u)j — u>o\,Xi <C ojq. The oscillating electric 
and magnetic dipoles of an isolated meta-atom radiate 
energy at respective rates Te and T 



M: 



r 



67reoC 3 
6nc 3 Lj 



(10) 
(11) 



resulting in the scattered fields E$ j and Hs j [see Eqs. (JSJ) 
and ([6])]. For simplicity, we assume that these radiative 
emission rates Te and Tm are independent of the res- 
onator j and that they are dominated by the meta-atom 
resonance frequencies, i.e. T-^j, Tmj <C CIq- We further 
assume the resonance frequencies occupy a narrow band- 
width around the central frequency of the incident field. 

The dynamics of current excitations in the meta-atom 
j may then be described by Qj(t) [introduced in Eq. ([I}] 
and its conjugate momentum 4>j (t) (with units of mag- 
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netic fluxi). In terms of the positive frequency compo- 
nents the equations of motion read^ 





*1 






r E 


)l 


— i — 

Uj 





-m, 



t) (12) 



E 



+ 

j,ext V 



(rj,t), (13) 



where the fields generated externally to meta-atom 
j that drive its dynamics, E+ ext (r, t) and H+ cxt (r,i) 
are produced by the sums of the corresponding inci- 
dent fields and the fields scattered by all other meta- 
atoms in the metamaterial sample, X^'=tj ^Sj'( T ^) anc ^ 
Ylj'^j Hs,j'( r ; i)i respectively. The component of the ex- 
ternal electric field Ej, cxt oriented along the dipole di- 
rection dj provides a net external electromotive force 
(EMF) hjdj • Et ext (rj,t) which drives cf>j(t). Similarly, 
the component of the external magnetic field Hj i0Xt along 
the magnetic dipole direction riij provides a net applied 
magnetic flux Aj^xhj ■ H~ h ext (r :) -, t) that drives Qj(t). In 
the absence of radiative emission and interactions with 
external fields, current and charge oscillate within the 
meta-atom at the resonance frequency ujj. The meta- 
atom dynamics are therefore naturally described by the 
slowly varying normal variables 



V2 



i(t) 



y/uj~Lj 



(14) 



In the absence of external field interactions and damp- 
ing, bj oscillates with frequency ((dj — i7o), i.e. bj(t) = 
b 3 (0) exp [— i (ujj — flo) t]. For nonzero Fe,Fm <C Qq, 
losses and driving from the external field act to perturb 
this oscillation. 

The current oscillation dynamics in the meta-atom j, 
described by Qj(t) and <f>j(t) in Eqs. (fl2|) and (fl~3|) . is 
driven by the incident field and the fields scattered from 
all the other meta-atoms and acts as a source of radi- 
ation that, in turn, drives the other meta-atoms. The 
expressions for the scattered fields by polarization and 
magnetization densities [Eqs. (O and ((6])] (generated by 
excitations in meta-atoms) and the expressions for the 
oscillating charge dynamics [Eqs. (fT2j) and (fT3| ] form a 
coupled set of equations, describing EM field mediated 
interactions between the resonators. In terms of the nor- 
mal variables bj, these interactions may be represented 
by the set of equations,— 

b = Cb + f in , 

where wc have defined 



(15) 



b(t) 



/ bi(t) \ 



\ b nN (t) J 



fin(t) 



/2,in(t) 



\ fnN,in(t) J 



(16) 



The driving fj t - m of each meta-atom j results from the 
EMF and magnetic flux induced by the incident fields^ 



The component of the incident electric field E; n par- 
allel to the electric dipole orientation dj induces the 
EMF, while the component of Hj n along the magnetic 
dipole orientation rhj provides an incident magnetic flux. 
Here we assume that the meta-atom magnetic dipoles are 
aligned perpendicular to the incident magnetic field, and 
thus only the EMF contributes to the driving of each 
meta-atom, which is given by 



=d i .E+(r i ,t). 



(17) 



The coupling matrix between the meta-atoms in Eq. (115 
reads 



C 



-iA 



r 



I + - (iC E + iC u + C x + Cj) 



(18) 



where I represents the identity matrix. Here the detun- 
ings of the incident field from the meta-atom resonances 
are contained in the diagonal matrix A with elements 



(19) 



and the energy carried away from individual meta-atoms 
by the scattered fields manifests itself in the decay rate 



r = r,= 



M 



r< 



(20) 



appearing in the diagonal elements of C. We account 
for non-radiative, e.g. ohmic losses, by introducing phc- 
nomcnological decay rate Tq. The multiple scattering 
processes are included in the terms Ce, Cm, and C x , 
which generate interaction between the meta-atom dy- 
namic variables The matrices Ce and Cm characterize the 
electric dipole-dipole and magnetic dipolc-dipole interac- 
tions, respectively. Additionally, the interaction embod- 
ied by C x arises from the electric field emitted by the 
magnetic dipole of one atom driving the electric dipoles 
of the others. Similarly, Cj results from the magnetic 
field produced by the meta-atoms' electric dipoles im- 
pinging on the magnetic dipoles of all the other meta- 
atoms. Because the interaction matrices Ce|m|x govern 
interactions between distinct meta-atoms, their diagonal 
elements are zero. Their off diagonal elements are given 
byi 

3 ~ 

-TEdj ■ G(rj -r r ,flo) -dj,, 



[Cu] jtj , 



r M thj ■ G(r 3 - , O ) ■ ihj 
tdj ■ G x (rj - rj,,n ) ■ rh^ 



(21) 
(22) 
(23) 



where T = \/TeJm is the geometric mean of the electric 
and magnetic dipole emission rates. 



III. ASYMMETRIC SPLIT RING RESONATORS 

In order to investigate the effects of inhomogeneous 
broadening of meta-atom resonance frequencies on a 



metamaterial's collective EM response, we consider an 
ensemble of asymmetric meta-molecules arranged in a 
regular lattice. To facilitate our description of this EM 
response, in this section we summarize the behaviour of 
a single ASR in the context of the model presented in 
Sec. ID 

An ASR is a variation on the split ring resonator used 
to produce bulk metamaterials with negative indices of 
refraction<22r— The meta-atoms of an ASR consist of two 
separate concentric circular arcs labeled by j G {1, r} 
and separated by u = r r — n. The current oscillations 
in meta-atoms produce electric dipoles with orientation 
d r = di = d associated with charge oscillating between 
the ends of the arcs. Owing to the curvature of the meta- 
atoms, these currents also produce magnetic dipoles with 
opposite orientations m r = — mi = rh where dlu and 
mlu,d. An asymmetry between the rings, e.g., result- 
ing from a difference in arc length, manifests itself as a 
difference in resonance frequencies with 



uj\ = ujq — Suj 



(24) 
(25) 



To analyze the dynamics of a single ASR unit-cell res- 
onator consisting of two meta-atoms, we apply the for- 
malism presented in Sec. [Til According to (TT5|) . the nor- 
mal variables b r and b\ that describe the current oscilla- 
tions in the right and left meta-atoms, respectively, are 
coupled by the EM fields so that they evolve according 
to 



_ r (ASK) 

bj J ~ c V h 



jr.': 
.//.J 



(26) 



Here £( ASR ) denotes the specific coupling matrix between 
the two meta-atoms that depends of the radiative elec- 
tric dipole - electric dipole, magnetic dipole - magnetic 
dipole, and electric dipole - magnetic dipole interaction 
processes between the two meta-atoms [See Eqs. I|18p- 
([23]) .]. On the other hand, the incident field produces 
the driving terms fj^ n for each meta-atom j = 1, r [Sec 
Eq. (HU).]. 

To analyze the modes of the ASR, we consider the 
dynamics of symmetric c+ and antisymmetric c_ modes 
of oscillation defined by 



1 



c \ 



71 (6r±6l) 



(27) 



The oscillations c± represent the eigenmodes of the ASR 
in the absence of asymmetry 5u> = 0. By diagonalizing 
C (ASR) with Su) = 0, one finds the eigenvalues of the 
modes c±, 



A± = -i (uj - O ± A) - -i 



(28) 



The interaction between the elements shifts the two 
collective resonance frequencies by equal and opposite 
amounts A and results in the decay rates j±, where the 
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FIG. 1. A schematic illustration of an array of ASR meta- 
molecules excited in the uniform phase magnetic mode. The 
meta-atom currents in each ASR oscillate with opposite 
phases, producing magnetic dipoles represented by the red 
arrows. 



coefficients A and j± depend on the radiative interac- 
tions between the two meta-atoms.—^ When the spacing 
between the arcs u <C A ( A = 27rc/f2o ), the decay rates 
simplify to 



7+ = 2r E - 
7- = 2r M 



r Q 
-r Q 



(29a) 
(29b) 



In this limit, the symmetric mode, possessing a net elec- 
tric dipole, emits electric dipole radiation, and the anti- 
symmetric mode, possessing a net magnetic dipole, emits 
magnetic dipole radiation. We therefore may refer to 
symmetric and antisymmetric oscillations as electric and 
magnetic dipole excitations, respectively 

A nonzero asymmetry, 8lo ^ 0, tends to couple the 
symmetric and anti-symmetric oscillations. One finds 
that, when driven by an external field, these oscillations 
in a single ASR evolve as^ 



c± 



-i (cj ± A - ) - ~y 



c± - i5ujc T + F± , (30) 



where the driving terms F± = (f 1% i n ±fi.i n )/^/2. The sym- 
metric and antisymmetric oscillations are driven purely 
by the electric and magnetic fields, respectively, and 
when the meta-atom separation F + oc d-E^R, t) 

and F_ cx m • B^(R, t), where R is the center of mass of 
the ASR. Therefore, an incident field with Ej n || d and 
Bj n -L iii only excites the symmetric mode when Slu = 0. 
However for Sco ^ 0, the asymmetry couples the symmet- 
ric and antisymmetric modes, and this incident field can 
resonantly pump the anti-symmetric magnetic mode via 
an effective two-photon transition^ 

In this article, we consider a 2D metamaterial com- 
prised of ASRs arranged in a regular array. The fields 
scattered from each ASR then mediate interactions be- 
tween them, resulting in collective modes of oscilla- 
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tion distributed over the array, each with its own res- 
onance frequency and decay rate. Figure Q] provides 
a schematic illustration of a mode consisting primarily 
of magnetic dipoles oscillating in phase throughout the 
metamaterial4£ We showed in Ref.ythat, for sufficiently 
large array, such a mode radiates more slowly than the 
magnetic excitation of a single ASR and is responsible 
for the transmission resonance observed by Fedotov et 
alX. When an incident EM field whose magnetic field is 
perpendicular to the ASR magnetic dipoles impinges on 
the array, this mode cannot be excited directly. But, the 
presence of an asymmetry provides a coupling between 
electric and magnetic dipoles allowing it to be driven. 
Driving of the metamaterial's collective modes are re- 
sponsible for the cooperative response that yields phe- 
nomena such as the observed transmission resonance*^ 



terms of the eigenvalue Aj by 



IV. INHOMOGENEOUS BROADENING 



In this section, we study the effects of inhomogeneous 
broadening on the cooperative EM response of an array 
of ASRs. Here inhomogeneous broadening refers to a 
statistical uncertainty in the resonance frequencies in in- 
dividual ASR meta-molecules. Such an uncertainty may 
result, for example, from imperfections in the manufac- 
turing processes which yield meta-atoms whose shape 
varies slightly from the design specifications. Collec- 
tive modes that are phase matched with an incident EM 
field have been shown to be responsible for transmission 
resonances^ that have been observed experimentally^^ 
We will illustrate the response of a regular array of ASRs 
to an incident plane wave, and show that broadening ad- 
versely affects the characteristics of the response respon- 
sible for observed resonances. 

We consider an ensemble of ASR meta-molecules 
whose constituent meta-atoms are separated by u = ue x 
with electric dipoles oriented along d = e y such that 
a symmetric oscillation in a single ASR produces an 
electric dipole along d, and an antisymmetric oscilla- 
tion produces a magnetic dipole along m = e 2 . The 
ASRs are arranged in an N x x N y 2D square lattice 
with lattice spacing a and lattice vectors ai = ae x and 
&2 = ae y . The sample is illuminated by a cw plane wave 
E^(r) = i£e a e lk ' r with k = ke z , coupling exclusively to 
the electric dipole moments of the ASRs. 

The vector of 2N x N y normal variables describing the 
state of current oscillations in each meta-atom obeys the 
coupled equations of motion in Eq. (|15j) . where the ma- 
trix C [See Eq. (|18|) .] arises from the meta-atom inter- 
actions mediated by the EM field. The metamaterial 
therefore exhibits 2N x N y collective modes of oscillation 
corresponding to the eigenvectors v ?; (i = 1 . . . 2N x N y ) of 
interaction matrix CA Each eigenmode i possesses a par- 
ticular resonance frequency fli and decay rate 7, given in 



Qi = -Im(Aj) + Q . 
7i = -2Re(Ai), 



(31a) 
(31b) 



respectively. 

Since the incident field drives all ASRs uniformly, it 
couples most effectively to the collective modes in which 
all of the metamolecules oscillate in phase. The two 
modes of particular interest are the uniform electric and 
uniform magnetic modes. In the absence of an asym- 
metry (5u> = 0), the incident field drives the uniform 
electric mode which, owing to the electric dipole ori- 
entations, emits strongly into the ±e z directions. The 
electric dipoles oscillating in phase are responsible for 
reflection from the metamaterial. By contrast, the mag- 
netic dipoles in the uniform magnetic mode, illustrated 
in Fig. [Tl emit into the plane of the meta-material array, 
and for sufficiently large lattices of closely spaced ASRs, 
the magnetic dipole radiation becomes trapped. This re- 
sults in a suppressed radiative decay rate of the uniform 
magnetic mode.— Introduction of an asymmetry (Scu 7^ 0) 
provides an effective coupling between these two collec- 
tive modes similar to the coupling between the symmetric 
and antisymmetric modes in a single ASR. The collective 
magnetic mode can thus be resonantly excited at the ex- 
pense of the electric dipoles resulting in a transmission 
resonanc o 5 ' 29 whose quality factor increases with the size 
of the array, as observed by Fedotov et aZ.— 

Specifically, the uniform magnetic mode is the eigen- 
mode v m of C [See Eq. (JTSJ) .] which maximizes the overlap 
O m (byi) with the pure magnetic excitation bA in which 
all meta-atoms arc excited with equal amplitude and the 
two meta-atoms in each split ring oscillate out of phase. 
Similarly, the uniform electric mode is the eigenvector 
v e of C that maximizes the overlap O c (bs) with the pure 
electric excitation bg for which all current oscillations os- 
cillate in phase with equal amplitude. Explicitly, these 
column vectors of 2N x N y elements are given by 



b A 



V 



-1 
1/ 



b s = 



1 



1 

Vi7 



(32) 



The alternating signs of the elements of 6a indicate the 
relative phase of the oscillations in each meat-atom of 



an ASR. We define the overlap of mode 
arbitrary excitation b as 



O m/e (b) 



Tm/c 



If 



EJvfbp 



with an 



(33) 



where the index i is summed over all the eigenmodes of 
the interaction matrix C. The uniform electric and uni- 
form magnetic modes for a 21 x 21 array of ASRs with 
a nonzero asymmetry parameter Slo = 0.3T are shown 
in Figs. [5] and [31 respectively We used the experimental 
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FIG. 2. The uniform electric mode v e of a homogeneously 
broadened ( a — ) 21 x 21 ASR square lattice, (a) The 
electric dipole excitations |c+,;| 2 and (b) the magnetic dipole 
excitations \c-j 2 of the uniform electric mode v e . The phases 
of the electric (c+ t i) and magnetic (c_,;) dipole excitations are 
indicated by the color of the surfaces in (a) and (b), respec- 
tively. The black dots indicate the positions of the ASRs in 
the array. The excitations | c-i- ,z | 2 are normalized to the peak 
ASR excitation Eg = max; (IC+,; | 2 + |c_j| 2 ). The vertical 
scale of panel (b) was amplified by a factor of ten to render 

the magnetic dipole excitations | c ,z | 2 visible. The vertices on 

the plots correspond to the ASR positions. The lattice spacing 
a — 0.28A, meta-atom separation within an ASR u — 0.12A, 
Te = Tm, and asymmetry parameter Suj = 0.3r. 

value for the lattice spacing from Ref. [H and the estimate 
the asymmetry parameter Suj ~ 0.3r from the relative 
arc lengths of the ASR meta-atoms studied by Fedotov 
et al.-^^iL The ohmic loss rate To was fitted so that the 
quality factor of the uniform magnetic mode as a function 
of system size matched the experimental observations £ 
Where the state of the ensemble is characterized by the 
vector of meta-atom normal variables b [See Eq. (jTSj)-]. 
the symmetric (electric) and antisymmetric (magnetic) 
oscillations of an ASR 1(1 = 1, ...N x N y ) are repre- 
sented by c +: i and c_j, respectively, where 

C±,, = -^(&2i-l±&2l) • (34) 

The respective symmetric (electric dipole) and antisym- 
metric (magnetic dipole) excitation energies in ASR I are 

proportional to |c+.;| 2 and | c ^ 1 2 - The asymmetry in 

the ASRs causes a mixing of the electric and magnetic 
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FIG. 3. The uniform magnetic mode v m of a homogeneously 
broadened ( a — ) 21 x 21 ASR square lattice, (a) The 
electric dipole excitations |c +i ;| 2 and (b) the magnetic dipole 
excitations |c-,j | 2 of the uniform electric mode v e . The phases 
of the electric (c+ ; j) and magnetic (c_,;) dipole excitations are 
indicated by the color of the surfaces in (a) and (b), respec- 
tively. The black dots indicate the positions of the ASRs in 
the array. The vertical scale of (a) was amplified by a factor 
of ten to render the electric dipole excitations c +i ;| 2 visible. 
All parameters are as in Fig. [2] 

dipolcs, producing a slight electric dipole excitation in 
the uniform magnetic mode v m and a slight magnetic 
excitation of the electric mode v e . In the example illus- 
trated in Figs. [Hand O when the ohmic loss rate in each 
meta-atom is To — 0.14r, the electric mode has an en- 
hanced decay rate j c = 2.7T and the magnetic mode has 
a suppressed decay rate 7 m = 0.3 IT with respect to the 
total isolated single meta-atom decay rate T. Due to the 
suppressed decay rate of v m and its small electric dipole 
component, this mode can be resonantly excited by the 
incident field. 

We illustrate this phenomenon in Fig. [51 where we show 
the steady state response [See Eq. (fl"5|).]. 

b r = -C-H in , (35) 

of an array whose resonance frequencies are not inhomo- 
geneously broadened. The right and left meta-atoms of 
ASR / in such an array have respective resonance frequen- 
cies u>2i-\ = Wo + &u an d 0J21 = ujq — 5uj centered around 
wo- The driving field is resonant on the uniform magnetic 
mode, and the asymmetry in the split rings Suj = 0.3T 
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FIG. 4. The response of a homogeneously broadened ( a = ) 
21 x 21 ASR square lattice to an incident plane wave elec- 
tric field resonant on the uniform collective magnetic mode 
v m , showing the excitation of a uniform magnetic dipole ex- 
citation at the expense of electric dipoles. (a) The electric 
dipole excitations |c+.;| 2 and (b) the magnetic dipole excita- 
tions |c_.;| 2 of the response. The phases of the electric (c +i ;) 
and magnetic (c—j) dipole excitations are indicated by the 
color of the surfaces in (a) and (b), respectively. The black 
dots indicate the positions of the ASRs in the array. The 
ohmic loss rate is To = 0.14T. All other parameters are as in 
Fig-H 



facilitates the phase-coherent excitation of the magnetic 
dipoles at the expense of the electric dipoles. Panels (a) 
and (c) of Fig. [4] illustrate that the magnetic dipoles are 
much more strongly excited than the electric dipoles in 
the bulk of the array, and that these magnetic dipoles 
oscillate in phase. The more excited of the weak elec- 
tric dipoles also oscillate in phase, thus facilitating the 
driving of this excitation by the uniform incident field. 
Although other collective modes of the system are ex- 
cited, more than 60% of the excitation energy resides in 
the uniform magnetic mode v m . In the absence of ohmic 
losses, one can optimize the asymmetry parameter Slu in 
large lattices so that over 98% of the excitation energy 
resides in the uniform magnetic mode£ 

The introduction of inhomogeneous broadening alters 
the collective interactions and can destroy the charac- 
teristics of the mctamaterial response that produces the 
transmission resonance. We model the inhomogeneous 
broadening by shifting the central resonance frequency of 
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FIG. 5. Response of inhomogeneously broadened square lat- 
tice to an incident plane wave resonant on the homogeneously 
broadened uniform magnetic mode v m . (a) The electric dipole 
excitations | c_t_ , z | 2 and (b) the magnetic dipole excitations 
|c_j| 2 of the response. The phases of the electric (c +i ;) and 
magnetic (c_,;) dipole excitations are indicated by the color 
of the surfaces in (a) and (b), respectively. The black dots in- 
dicate the positions of the ASRs in the array. The individual 
ASR resonance frequencies are shifted by a independent iden- 
tically distributed Gaussian random variables with standard 
deviation a — O.SSuj. All other parameters are as in Fig. [4] 



each ASR I by independent identically distributed Gaus- 
sian random variables X\ with zero mean and standard 
deviation a. With the asymmetry characterized by <5w, 
the right and left circular arcs in ASR I possess resonance 
frequencies LO21-1 = luq + Xi + Slo and LO21 = + Xi — Slu, 
respectively. The deleterious effects of this broadening 
on the response are illustrated in Fig. [5l which shows a 
much less uniform magnetic response in addition to lo- 
calized electric dipole excitations. This non-uniformity 
inhibits the coherent reflection and transmission through 
the meta-material array. 

We quantify the effects of inhomogeneous broadening 
by examining the overlap of the metamaterial steady 
state response b r with the magnetic mode v m (er = 0) of 
a homogeneously broadened array (corresponding to the 
case in which all ASR meta- molecules are identical). If 
the excitation is purely in the mode v m (0), then the over- 
lap O m (b r ) [Eq. (|33l)] is unity. Figure^a) shows the over- 
lap O m (b r ) of the response to an incident field resonant 
on the homogeneously broadened mode v m (cr = 0) aver- 
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a/Sco 5 a/8co 5 

FIG. 6. (a) The overlap of the uniform magnetic mode with 
the excitation driven by an incident EM plane wave and (b) 
the decay rate 7m(c) of the uniform magnetic mode v m (a) 
as a function of inhomogeneous broadening a. These quanti- 
ties are plotted for Suj = 0.3r and ohmic losses To = 0.14F 
(solid black line), Suj = 0.3F and Fo = (dashed blue line), 
and Suj = 0.1F and To = (dot dashed-red line). The lines 
indicate the average over 240 sample realizations, and the er- 
ror bars indicate the standard deviations. The incident wave 
has electric field polarization e y aligned with the ASR elec- 
tric dipoles and is resonant on the mode v m (a = 0) of the 
homogeneously broadened array. The decay rates of the uni- 
form magnetic modes for a homogeneously broadened array 
7m(o- = 0) are 0.308F (black line), 0.145F (blue line) and 
0.037ir (red line). This shows that excitation of the mag- 
netic mode, and hence the transmission resonance, vanishes 
as the inhomogeneous broadening becomes comparable to Suj. 



netic mode more difficult to excite since any excitation of 
this mode is more quickly radiated away. In order to for 
an array of ASRs to exhibit a transmission resonance, a 
large fraction of the excitation created by the driving field 
must be in the uniform magnetic mode distributed over 
the array. In low loss metamaterials, this can be achieved 
in conjunction with a higher quality for that resonance 
for larger arrays and smaller values of Suj However, as 
Fig. [6] indicates, reduction in the asymmetry to achieve 
this quality factor enhancement correspondingly reduces 
the tolerance for inhomogeneous broadening in the reso- 
nance frequency. 

The observation that the cooperative metamatcrial re- 
sponse to EM fields can be suppressed in the presence 
of sufficiently strong inhomogeneous broadening could 
potentially also be exploited in design of metamaterial 
samples that would benefit from well-defined homoge- 
neous properties for electric susceptibility and magnetic 
permeability, such as diffraction-free lenses due to nega- 
tive refractive indexi^^r— One could prepare a controlled 
amount of inhomogeneous broadening for the metama- 
terial sample in order to generate an EM response that 
more closely mimics standard continuous medium elec- 
trodynamics with suppressed contribution from recurrent 
scattering events. 



aged over 240 realizations. The solid black line was calcu- 
lated for the same parameters as in Figs.[2]through[5]with 
varying degrees of broadening. The blue dashed curve 
shows the corresponding overlap in the absence of ohmic 
losses, while for the red dashed-dot curve, the asymmetry 
parameter was reduced to Suj = 0.1 and Tq = 0. In all 
cases, as the broadening standard deviation a becomes 
comparable to Sui, the ability to excite the uniform mag- 
netic mode drastically decreases. The other modes that 
are excited either contain electric dipole components or 
are not phase matched. This either results in scattering 
of the field or in absorption of the field due to ohmic 
losses. As a result, the coherent collective response re- 
sponsible for the transmission resonance observed in by 
Fedotov et a& becomes unobservable when the inhomo- 
geneous broadening is larger than Suj. These collective 
effects do persist, however, for a roughly half Suj. Wc 
show the effect of inhomogeneous broadening in the de- 
cay rate 7 m (c) of the magnetic mode itself in Fig. [6^b) 
as a function of broadening. Randomization of the ASR 
resonance frequencies apparently has little effect on the 
collective lincwidth of the magnetic mode for a < Suj. 
However, for larger degrees of broadening, the decay rate 
of this mode can be increased several times over and the 
resonance lincwidth narrowing, that results from the co- 
operative response of the metamaterial array, disappears. 
Furthermore, the large standard deviations of 7 m (c) in- 
dicate that the width of the uniform magnetic mode is 
highly sensitive to the particular realization of ASR res- 
onance frequencies. A larger decay rate renders the mag- 



V. CONCLUSION 

In conclusion, we have analyzed the collective modes of 
a finite-sized 2D metamaterial array of ASR resonators 
and how they are influenced by an inhomogeneous broad- 
ening of the resonance frequencies of the individual res- 
onators. The study was motivated by recent experimen- 
tal observations of transmission resonance linewidth nar- 
rowing as a function of the size of the system.— This 
effect can be understood by analyzing the resonance 
linewidths of collective modes of the system that undergo 
dramatic narrowing due to strong EM field mediated in- 
teractions between the resonators. As demonstrated by 
previous comparisons between the numerical simulation 
results and the experimental observations,— the response 
can be analyzed by a simplified model in which each 
meta-atom is treated as a discrete element supporting 
a single mode of current oscillation possessing electric 
and magnetic dipole moments. Collective interactions 
between the meta-atoms are mediated by the scattered 
EM fields. The excellent agreement between the theory 
and the experiment can be understood by a relatively 
weak higher-order multiple radiation of individual ASR 
metamolecules .— 

We examined in detail how inhomogeneous broadening 
of resonator resonance frequencies impairs the coherent 
collective phenomena that are expected to find impor- 
tant applications in metamaterial systems i & ' 9 ' 29 While 
the transmission resonance experimentally observed in 
Ref. [l| persists for inhomogeneous broadening that is a 
fraction of the ASR asymmetry parameter Suj, the co- 
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operative response vanishes when the broadening begins 
to exceed that parameter. Production of high quality 
resonances with low loss materials requires the reduction 
of du!& Figure [H] indicates that, in order to produce such 
high quality resonances, the uniformity in the production 
of meta-molecules will need to become correspondingly 



small. 
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